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The possibility of determining the state of a quantum sys- 
tem after a continuous measurement of position is discussed 
in the framework of quantum trajectory theory. Initial lack 
of knowledge of the system and external noises are accounted 
for by considering the evolution of conditioned density matri- 
ces under a stochastic master equation. It is shown that after 
a finite time the state of the system is a pure state and can 
be inferred from the measurement record alone. The relation 
to emerging possibilities for the continuous experimental ob- 
servation of single quanta, as for example in cavity quantum 
electrodynamics, is discussed. 

42.50.Lc,03.65.Bz,42.50.Ct 



I. INTRODUCTION 

There is currently great interest in experiments which 
obtain useful information about a quantum system or 
state in single runs of an experiment. Very recently it 
has been possible to distinguish the quantum mechanical 
and classical models of the interaction of an atom with a 
mode of a high finesse optical cavity as the result of con- 
tinuous monitoring of the output light while a single atom 
passes through the cavity (lj . As a result of this contin- 
uous monitoring the quantum mechanical backaction of 
the measurement process may be expected to have a sig- 
nificant effect on the evolution of individual runs of the 
experiment. Moreover, it may be possible in the future 
to modify the evolution of the system through feedback 
based on this continuous measurement Current ex- 
perimental technology, such as that described in Q], is 
reaching the point at which determining the state of the 
system and observing the effects of backaction and feed- 
back in a single run of an experiment is a real possibility. 
With this situation in mind, we discuss the identification 
of the state of the system following a period of contin- 
uous observation and the extent to which this state can 
be tracked, taking into account factors such as imperfect 
initial knowledge of the state and imperfect detection ef- 
ficiency. We focus on a model of continuous position 
measurement of a mechanical oscillator which is relevant 
to the experiment of Mabuchi et al. but which also has 
relevance to other instances of interferometric position 
monitoring such as gravitational wave detection. 

The problem of describing quantum systems undergo- 
ing continuous measurement has attracted much theoret- 
ical interest in recent years. As discussed by Wiseman 
P| these theories admit a variety of interpretations; as 



tools for efficient stochastic calculation of ensemble aver- 
ages in lieu of solving master equations [Q , as equations 
describing the evolution of systems conditioned on mea- 
surements and as a description of the evolution of 
a system coupled to an environment, in which collapse 
of the wavefunction is supposed to be associated with 
the coupling to the environment ||. Here, we take the 
second viewpoint namely that the conditioned state rep- 
resents the observer's best description of the system state 
given the results of the continuous measurement process. 
Adopting the first or third viewpoints one is led to de- 
scribe the system by a pure state vector throughout the 
evolution although the reasons for doing so are some- 
what different in each case. By contrast, a description of 
one's conditioned state of knowledge necessarily requires 
mixed states in order to account for incomplete knowl- 
edge of the system. From this viewpoint the fundamental 
equation for the conditioned evolution is the stochastic 
master equation (SME) This is able to account for 
the effects of mixed initial states, imperfect detection effi- 
ciencies and the existence of unmeasured couplings to the 
environment. However, to date, relatively little work has 
attempted to address the evolution of the conditioned 
state in any of these situations po| . In this paper we 
consider a system which is simple enough that almost 
all the work can be done analytically and which admits 
a treatment of all of these imperfections. This helps in 
developing intuition about the role of SME's and their 
possible relevance to experiments. 

A projective measurement has the property that if the 
result of the measurement is known, the state after the 
measurement is pure, and depends only on the measure- 
ment result. It would be hoped that in a continuous 
measurement there would be some finite interval of time 
after which the measurement has effectively given rise to 
a projection, so that the system is placed in a particular 
state which depends only on the sequence of measure- 
ment results and which can be calculated without knowl- 
edge of the initial state. If the resulting state is pure 
then a stochastic Schrodinger equation (SSE) would be 
a perfectly adequate tool for describing the subsequent 
system evolution. In this paper we investigate the con- 
ditions which lead to such an effective collapse and over 
what timescale it takes place. This is made possible by 
considering density matrices and the SME rather tha 
wave vectors and the SSE. In a real experiment there 
will also be uncontrolled, unmeasured couplings of the 
system to the environment and in this case the effects 
of the measurement will compete not only with the co- 
herent internal dynamics of the system but also with the 



1 



randomizing effects of the coupling to the bath. This 
may lead to mixed conditioned states even after long pe- 
riods of continuous measurement and limit the observer's 
ability to make inferences about the system state. Un- 
derstanding the process by which the conditioned state 
may collapse onto a pure state and the effects of noise as 
described by the SME allows us to define conditions un- 
der which continuous measurements in real experiments 
are approximations to ideal measurements. 

This paper is organized as follows. In Sec. [n] we estab- 
lish our simplified model of the continuous position mea- 
surement of an oscillator and solve the SME for Gaussian 
initial states. We find the time over which the second- 
order moments approach their steady-state values and 
calculate the entro py of the conditioned state as it be- 
comes pure. In Sec. Ill we discuss the classical problem of 



state identification for the noisy measurement of the the 
position of an oscillator and derive a kind of uncertainty 
principle relating the observation and process noises if 
the classical model is to reproduce the SME. In Sec. IV 
we show that the time-scale over which the second-order 
moments of the conditioned state reach their steady state 
is the same as that over which the conditioned state is 
completely determined by the measurement record. Sec- 
tion ^ discusses the effect of heating, noise and detection 
inefficiency on these conclusions. Finally in Sec. ^ we 
summarize and make some comments about future ex- 
tensions of this work. 



II. SOLVING THE STOCHASTIC MASTER 
EQUATION 

A. A Generalized Model for Continous Position 
Measurement 

In this paper we shall consider the abstract model of 
continuous position measurement discussed by Caves and 
Milburn 11 |] with a harmonically bound rather than 
a free measured particle. Projective position measure- 
ments are imagined to be made on a sequence of meters 
coupled briefly to the system and the limit of very fre- 
quent meter interactions and very broad initial meter po- 
sition distributions is taken. This leads to a continuous 
evolution for the system of interest. Although this model 
should correspond in some limit to any continuous posi- 
tion measurement of a single oscillator at the standard 
quantum limit, one system which does realize it at least 
approximately is the dispersive regime of single atom cav- 
ity quantum electrodynamics (QED). In this system the 
position of an atom inside a high-finesse optical cavity 
causes a phase shift of the field driving the cavity which 
can be monitored by homodyne detection |l^-|l5|]. The 
SME for the conditioned evolution of the system in its 
Ito form is [IlT 



dp c = -i[H sys , p c ]dt + 2aV[x]p c dt + V2aH [e'^x] p c dW. 

(2.1) 

The superoperators T>[c] and 7i[c] acting on a density 
matrix p are T>[c]p — cpc< — ^c^cp — ^pc'c and Tt[c]p = 
cp + pc) — Tr(cp + pc*)p where c is an arbitary operator. 
We will imagine that the atom is harmonically bound, 
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(2.2) 



The constant a describes the strength of the measure- 
ment interaction and in the cavity QED example depends 
on the strength of the coherent driving and of the damp- 
ing of the cavity. For the moment we will consider a one- 
sided cavity and perfect detection so that all the output 
light is detected. This assumption will be relaxed in Sec. 
(|v|) . In a generalization of the Caves and Milburn model 
we will allow projective measurements of any quadrature 
of the meters, not just position, since this can be realized 
in the cavity QED experiments by varying the phase <p 
of the local oscillator in the homodyne detection of the 
output light. The resulting measurement current i = 5? 
(suitably scaled) is 



dQ = cos((/))(x) c dt + \ —dW. 

8a 



(2.3) 



This stochastic master equation with the full depen- 
dence on <fi was discussed by Diosi |l7j in the context 
of a phenomenological model of position measurement 
through photon scattering where the kind of measure- 
ment made on the scattered photon determines the value 
of (p. Clearly if we choose (j> = 0, ir the homodyne detec- 
tion is an effective measurement of the atomic position, 
whereas if <j) = 7r/2,37r/2 the measurement results are 
independent of the system state and only contain infor- 
mation about the noisy potential seen by the atom. For 
<f> = 7r/2,37r/2 the conditioned evolution is linear in the 
system state [^,|l5[. This is somewhat like a continuous 
quantum eraser where the continuous measurement of 
one quadrature of the probe system destroys the posi- 
tion information written into the other quadrature of the 
probe. For cf> = the SSE corresponding to this model in 
the case of pure states of the system has been considered 
in p^Jl9[ |. In this work we will use a simpler means of 
solving the evolution equations which straightforwardly 
applies to mixed states. 

It has been shown by Jac obs and Knight jl9| that the 
SSE corresponding to Eq. (2T) is one for which Gaus- 
sian pure states p0| remain Gaussian and pure under 
the evolution. Thus, if the system is initially in a mix- 
ture of Gaussian pure states, the c ond itioned state will 
remain Gaussian under the SME (2.1). This property 
holds true for single mode systems where the Hamilto- 
nian is at most quadratic in x and p and the operator c 
appearing in the Lindblad term T>[c] is linear in x and p 
and in all likelihood for multimode linear systems also. 
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- (x),(p),((Ax) 2 ) = (x 2 ) - (*) 2 ,((Ap) 2 > 
and (AxAp) S y lrL = \{{xp + px) - 2(x)(p)). 



If we restrict ourselves to Gaussian initial states — for 
example to thermal states of the oscillator — then there 
are only five quantities which completely define the state 

= (P 2 ) - (P) 2 
From now on 

we will use (c) to indicate the conditioned expectation 
value Tr(cp c ). The requirement that the initial states be 
Gaussian is not unduly restrictive since these are typi- 
cally the states that are most stable in linear systems 
and would therefore be a reasonable assumption for an 
initial state. Moreover, there is numerical evidence that 
non-Gaussian pure initial states become approximately 
Gaussian under the stochastic Sch rod inger equation evo- 
lution corresponding to the SME (2.1) on timescales fast 
compared to those considered here | 2l| , p2| . It has also 
been shown that in at least one such linear system an ar- 
bitary density matrix can eventually be written as a prob- 
abilistic sum over Gaussian pure states after sufficient 
evolution subject to the unconditioned master equation 
p3j . In considering such a linear system we are in effect 
specifying a semiclassical evolution, since the equations 
of motion for the Wigner function of the state in the un- 
conditioned evolution are similar to the classical Liouvillc 
equations for a phase space density, see for example [^4j . 
However we here include the important quantum feature 
of the measurement backaction, which is represented by 
the fact that the momentum diffusion is determined by 
the accuracy with which the particle's position is moni- 
tored So, as the noise in the position measurement is 
decreased by increasing a, the momentum diffusion ex- 



pressed by the Lindblad term in the SME (2.1) increases. 

Following Breslin and Milburn (2^] we can de- 
rive a system of differential equations for th e first 
and second-order moments from the SME (2.1) since 
d(c) =Tr(cdp c ). A similar calculation is performed 
in [p3[ . We define the dimensionless quantities x — 
(x) I y // h/2muj, p — (p) J ^hmuj/2, the second-order mo- 
ments V xx = 2muj((Ax) 2 )/h, V pp = 2((Ap) 2 } /hmuj, 
V xp = 2(AxAp) synl /h and a dimensionless parameter de- 
scribing the relative strengths of the measurement and 
harmonic dynamics, r = muj 2 /2ha. The Heisenberg un- 
certainty principle now requires that V XX V PP > 1 



state has the property that V XX V P} 



> 1. A pure 
xx ' P p - V 2 p = 1, represent- 
ing the fact that a Gaussian pure state is a minimum 
uncertaintly state for some pair of conjugate quadrature 
variables ^0|. In terms of these new quantities, the Ito 
stochastic differential equations for the first and second- 
order moments and the measured photocurrent are 



(2.4a) 



dp = -ujxdt + \ — (cos((f)V xp - sm(4>)) dW, (2.4b) 



dx = ujpdt + \ — cos(<fi)V xx dW, 



dQ = cos (/ft) xdt -, . 

V 2w 

1 dV XX 2 2/ ; , T -■ 

"77 = 2V xp ~ - COS Z ((j))V a 

uj dt r 



dW, 



2 

xx ' 



(2.4c) 
(2.4d) 



1 dVpp „/ sm(2(f>)\ Tr 2 2/I . 

= -2 1 — )V xp + - cos 2 </> 

uj dt \ r J r 

~^cos 2 (0)^ 2 p , 



ldK 
uj dt 



2 

--cos (4>)V xx V xp . 



(2.4e) 



(2.4f) 



As in |2q| , the Ito rules for stochastic differential equa- 
tions ]26|] and the properties of Gaussian states (2^] result 
in deterministic equations for the conditioned second- 
order moments which are decoupled from the equations 
for the means. The constant term in the equation for 
V pp refers to the momentum diffusion due to the position 
measurement and remains in the master equation for the 
unconditioned evolution. The non-linear terms describe 
the conditioning of the state on the measurement. The 
noisy contribution to the equation for dx seems a little 
like a stochastic impulsive force, however it is perhaps 
better to think of this term as updating the expected po- 
sition given the measurement result dQ in analogy with 
classical Bayesian state estimation. 

Equations like those above for the second-order mo- 
ments of the conditioned state arise very frequently in 
classical continuous-time observation and control prob- 
lems. They can be collected into a Riccati matrix differ- 
ential equation [p8| for the covariance matrix, 



— V = uj (C - VBV - DV - VA) , 
dt 



V = 



Vxx ^xp 

Vxp Vpp 



i(20) 



A = D 1 = 



_ [ 1 _ sin ( 2< ft) ) o 



B 

C 



2 cos (j) q 
r 

o o 
o o 

q 2 cos 



(2.5a) 
(2.5b) 

, (2.5c) 

(2.5d) 
(2.5e) 



A single variable Riccati equation which arose from the 
stochastic Schrodinger equation for this system was found 
and solved in |^8| . 

In practice it may not be homodyne but rather het- 
erodyne detection which can be experimentally achieved 
with noise at the quantum limit Q. In this case, the 
local oscillator is detuned from the cavity frequency 
by a frequency Ahct which is large compared to all 
system frequencies, with the result that the phase <j> 
changes very rapidly. The quantum theory of hetero- 
dyne detection is described by Wiseman and Milburn 
fl29| . The appropriate conditioned evolution can be de- 
scribed by averaging all the tr igon ometric functions of 
cf> in the evolution equations (2.4) except where they 
are multiplied by Ito increments. Thus the equations 
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for the second order moments are exactly those for ho- 
modyne detection with = where r is replaced by 
2r, corresponding to halving the signal to noise ratio 
of the measurement. Considering the stochastic inte- 
grals J t ' +5t cos(A hct t')dW(t')J* +St sm(A hct t')dW(t'), in 
the limit firstly of infinite Ahct then of infinitesimal St, 
leads to equations for the first order moments under het- 
erodyne detection 



dx = Lopdt - 



V xx dW u 



(2.6a) 



dp = -uxdt + J - (V xv dW x - dW 2 ) , (2.6b) 

where dW\ and dW 2 are independent Wiener increments. 
Again this is formally identical to the evolution with ho- 
modyne detection and = where r is replaced by 2r 
and in which a second independent noise process dW% is 
present. Scaled versions of the two quadrature compo- 
nents of the experimental photocurrent are given by 



(2.7a) 
(2.7b) 



dQi = xdt + A —dWi, 

V LU 

dQ 2 = dW 2 . 



Again note the replacement of r by 2r in the equation 
for dQi as compared to the equation for dQ with = 0. 
So if the quadrature-phase I 2 (t) current is collected and 
used to account for the noisy potential, or alternatively 
fed back in order to compensate this evolution, then het- 
erodyne detection is equivalent to homodyne detection 
with half the signal to noise ratio as far as the motional 
state is concerned. 



B. Steady State Conditioned Variances 

For all phases of the local oscillator ^ 7r/2,37r/2 
the second-order moments possess a steady state. For 
example if = 




(2.8a) 

(2.8b) 
(2.8c) 



which defines a pure state that agrees with the solution 
given in |is|| . The steady states are found to have ex- 
actly the same second order moments regardless of the 
initial purity of the system. Assuming ideal detection, 
the observer is eventually able to ascribe a pure state 
to the system. When the harmonic oscillator dynam- 
ics dominates over the measurement (r 3> l)the steady 



conditioned state is approximately a coherent state with 
Kra — ^pp — 1' Vxp — l/ r - K the measurement dynamics 
dominate (r < 1) then V™ ~ ^/r, V™ ~ V™ ~ 1 and 
the conditioned state is strongly squeezed in position, as 
one would expect for a position measurement which is 
rapid enough to overcome the internal dynamics of the 
system. The product of the position and momentum vari- 
ances is greater than that required by the Heisenberg un- 
certainty principle as a result of the internal Hamiltonian 
which gives a non-zero correlation V xp . Finally, the scal- 
ing we have chosen for the variances makes the limit of 
the free particle appear singular, however this limit exists 
and the results agree with Belavkin's |fisj| . Pure steady 
states also exist for other values of but the full ex- 
pressions are rather complicated so we will just consider 
two special cases. If the oscillator dynamics dominate 
(r> 1) the steady conditioned states are insensitive to 
the local oscillator phase leaving the conditioned state 
nearly in a coherent state V* x — 1 + sin(20)/2r, ~ 
1 — sin(20)/2r, ~ cos 2 0/r. If r <C 1 the conditioned 
state is strongly dependent on the choice of local oscilla- 
tor phase, 



■yl/ |cos0| + tan0, 



V™ ~ -^VV|cos0| +tan0, 



V*p ~ 1/ |cos0| +tan0. 



(2.9a) 

(2.9b) 
(2.9c) 



Rather surprisingly it is possible for some phases of the 
local oscillator that the momentum variance is in fact 
smaller than the position variance. This is because the 
measurement for non-trivial is a simultaneous measure- 
ment of the position of the oscillator and the momentum 
kicks to which it is being subjected and this can result in 
a more sharply defined momentum than position. This 
is really only a possibility for phases of the local oscilla- 
tor where there is very little position information in the 
record. 

In the cases = 7r/2,37r/2 the differential equations for 
the second order moments of the state are simply those 
that result from the unitary evol ution of a simple har- 
monic oscillator and the SME (2.1) describes a stochastic 
unitary evolution. As a result there is no steady state for 
the moments, which grow according to the conventional 
Schrodinger equation. In this case the measurement cur- 
rent is white noise. In |H| it was noted that where a 
Lindblad operator is hermitian there exists an unravel- 
ling of the master equation which does not localize the 
conditioned state. The reason for this is clear in this con- 
text, such an unravelling corresponds to a measurement 
in which the observer obtains no information about the 
system state. 
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C. Timescale for Determination of a Pure State 
through Measurement 



The matrix Riccati equation (2.5a) has an analytic so- 
lution given by Reid p8[ . Where U(t) and Wit) obey 
the linear coupled matrix equations 



— U = AU + BW, 
dt 

4-W = CU -DW 
dt 



(2.10a) 
(2.10b) 



and for times where U(t) is non-singular, the solution 
for the covariance matrix is V(t) = W(t)J7 _1 (t). The 
full solution for all the second-order moments and ar- 
bitary initial conditions is complicated and not particu- 
larly illuminating. In order to expose the general form 
of the solution we will just consider the position vari- 
ance in the case of an initial state of the oscillator with 
V xx (0) = V pp (0) = Vo, V xp (0) = 0, 



V xx (t) = 



V (c 2 cosh 2bt + b 2 cos 2ct) 
(1 + V 2 ) (c sinh 2bt - b sin 2ct) 



(b 2 + c 2 




V (c a sinh 2bt 
sinh 2 bt 



b 6 sin 2ct) 
b 2 sin 2 ct) 



1/V SS 



(2.11a) 



(2.11b) 



(2.11c) 



We have scaled time by the harmonic oscillation fre- 
quency u>. When 2bt ^> 1 the system is close to its 
steady state value regardless of Vo. The non-linearity 
of the terms describing the conditioning of the system 
state cause the time over which the conditioned state be- 
comes pure to be independent of the initial temperature 
of the state. For dehniteness we will define a collapse 
time T C oi = 2/buj as being the time at which the state has 
become effectively pure. When r > 1 this collapse time 
is r co i ~ 2r/u! — muj/fia since in this regime b ~ 1/r. In 
this regime c ~ 1 + l/2r 2 and there are many oscillations 
of the particle before the state is determined. That the 
time to determine a particular pure state of the system 
should increase with the frequency seems reasonable since 
unexpected values of the measurement current could be 
due to a mistaken idea of the position of the particle, 
the white noise in the measurement record or to motion 
due to the oscillation and these possibilities will be more 
difficult to distinguish if the atomic motion is fast. For 
smaller r, the measurement is becoming very good and 
this estimate for the collapse time is optimistic since it 
is hard to determine the state of the system in less than 
one period of the mechanical oscillation. Reductions in 
the conditioned momentum variance will only occur as 
the Hamiltonian evolution creates a correlation between 
the position and the momentum of the state. Even if the 



harmonic potential were absent a continuous measure- 
ment of position will give some information about the 
momentum of the particle. When r <C 1 the particle is 
essentially free as far as the measurement is concerned 
and the time for the state reduction to occur turns out 
to be r co i ~ ^J&m/fia which is determined solely by the 
strength of the measurement and the mass of the particle. 
In this situation b ~ c ~ l/y/r. Increasing the measure- 
ment coupling means that the time for the measurement 
to take place is reduced and in the limit of infinite a the 
model essentially describes a projective measurement of 
the position. If heterodyne detection is used rather than 
homodyne detection then r is replaced by 2r in the above 
equations with obvious implications for the timescale of 
the system collapse. 

A pure state describes a situation in which an observer 
has maximal information about the system. A mixed 
state describes less than maximal knowledge of system 
and it turns out that the amount of missing information, 
in an information theoretic sense, required to complete 
the specification of the state may be measured by the von 
Neuman entropy S(p) = — Tr(plogp) of the density op- 
erator, see for example |53|. Thus the entropy allows us 
to quantify the extent to which the measurement has de- 
termined the state of the system at a given time and also 
the extent to which other environmental couplings limit 
what an experimenter can say about the system state. 
Another commonly used measure of the "mixedness" of 
a given density matrix is the linear entropy or purity 
s(p) = 1— Tr(p 2 ). For a single mode Gaussian state these 
quantities are simple functions of the unitless "area" of 

the state in phase space A - 



V XX V PP - V 2 p |31 



1 

A' 



(2.12a) 



S(p) = £±i m (A + l) - £—L\n(A- 1) - In 2. (2.12b) 

Note that A is just the determinant of the covariance 
matrix of the position and momentum probability dis- 
tributions for the conditioned state. For a pure state 
A = l,s(p) — > 0, Sip) — > and as the state becomes 
increasingly mixed it occupies a larger phase space area 
such that as A — > oo, sip) — > 1, S(p) — > oo. As we would 
expect if the state is widely spread in phase space then 
our knowledge of the system is poor and the information 
needed to complete the description of the state is large. 

The time evolution of the variances and the linear and 
von Neumann entropies of the conditioned state with 
14.(0) = VpJO) = 20,V xp = 0,u = l,r = 20 is plot- 
ted in Fig. (|l]). These parameters are chosen since the 
measurement dynamics are not fast enough to obscure 
the harmonic oscillation totally and because achieving 
very small values of r will probably be difficult in prac- 
tice. Several features of Fig. (jl]) are relevant. Firstly the 
initial very rapid reduction of the position variance is as- 
sociated with the the first part of the measurement record 
making a reasonably accurate determination of the posi- 
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tion. Then over the timescale of the harmonic oscillation 
the momentum variance also reduces as the dynamics 
correlate the position and momentum. Note that the re- 
duction in momentum variance occurs only when there 
is a strong correlation between the position and the mo- 
mentum. As V xp becomes small the position variance 
decays more rapidly and the reduction of the momentum 
variance slows. Eventually all the second-order moments 
decay to the steady state values predicted above. This 
initial fast reduction of the position variance is accom- 
panied by a fast reduction of the von Neumann entropy 
which damps to zero as the system approaches steady 
state. 



D. Cavity QED Realization 

Although we have been considering an abstract model 
for continuous position measurement this work is moti- 
vated by emerging experimental possibilities in areas such 
as cavity quantum electrodynamics (QED). The position 
dependent phase shift induced by an atom strongly cou- 
pled to a high finesse optical cavity mode in the dispersive 
limit of cavity QED [[l2|-|l5| realizes the abstract position 
measurement coupling considered here. It is currently 
possible to detect the presence of a single cold atom in the 
cavity through measurements of the output field [ f32|]33| 
and great progress has been made in observing single 
atom events with broad bandwidth close to the disper- 
sive regime The phase shift changes most quickly 
with position where the gradient of the field is greatest 
and if the atom is harmonically confined in this region 
then the model discussed above would be approximately 
realized. There is also the far off-resonant optical poten- 
tial which will lead to large forces on the atom in this 
regime which would move the atom quickly away from 
this region of the standing wave. However an optical 
standing wave at a nearby frequency could in principle 
be tuned to cancel the ac-Stark shift of the ground state 
in the region of the harmonic confinement. In fact the 
dipole force could straightforwardly be included in our 
simplified model resulting in only minor changes for suf- 
ficiently strong harmonic confinement. For the moment 
we will not specify the exact source of the potential con- 
fining the atom, but the use of light forces from a far 
off-resonant optical standing wave or a standing wave in 
another mode of the cavity are possibilities. There is also 
experimental work aimed at confining ions in high finesse 
optical cavities which could realize such a system p3 . 

If we imagine a harmonic potential confining the atom 
to this region of the standing wave with some constant 
restoring force to overcome the dipole force then in the 
far-detuned and Lamb- Dick e limit the resulting SME 
would be exactly Eq. (2.1) above pj|. The constant 



a would then be equal to 2g(pk L /A k, where go is the 
maximal single photon Rabi coupling in the cavity, n is 
the number of photons present in the driven cavity in 



the steady state, A is the detuning between the atomic 
and cavity resonances — the external laser driving is on 
the cavity resonance, n is the cavity field decay rate 
and ItL is the wave number for the light resonant in- 
side the cavity. Using the parameters of Jl5| which are 
based on cavity parameters achieved by Hood et al. |33j | 
gives a = 2.4 x 10 20 s _1 m -2 and this determines the 
rate of decay of the off-diagonal terms in the position 
representation of the density matrix under the uncondi- 
tioned evolution. This rather large number means that 
the density matrix elements (^Ipja;') where x and x' are 
separated by nine nanometers or around 1% of a wave- 
length will decay in the unconditioned evolution at the 
rate 2.0 x 10 4 Hz. The decay of off-diagonal elements 
of the density matrix in a particular basis is often asso- 
ciated with decoherence and the emergence of classical- 
ity In this case the decoherence is due in to the 

measurement coupling. In jl3| this decay of off-diagonal 
density matrix elements is described as state reduction, 
in this work we are interested in state reduction onto a 
pure state and the rate of this process is determined not 
only by a but also by the length scale of a typical pure 
state of the uncoupled system. Thus we found above 
that reduction onto a pure state took place in a time 
t co i ~ mu/Aha and the dependence on the length scale 
of the harmonic oscillator is clear. In order to find the 
rate of collapse onto a pure state in the conditioned evo- 
lution it is therefore also necessary to know the oscillation 
frequency of the atom due to its harmonic confinement. 
Assuming this is achieved optically, the value for lu/2tt 
could be in the range of tens to hundreds of kilohertz. 
So for example in |15| the potential for a cesium atom 
resulting from the same cavity and driving parameters 
gives uj/2-k = 2 y /hkl/2my / g : n/A/2Tr = 180 kHz, while 
lo/2tt ~ 60 kHz has been achieved for cesium in optical 
lattices by Haycock et al. |3(| . For such a hypothetical ex- 
periment with cesium we now have r = 5.6 and r = 0.63 
respectively. Estimates, as outlined above, for the time 
for an experimenter to determine a pure state of the 
system through heterodyne detection are then 19ps and 
8.9/iS. Both these times are reasonably close to the min- 
imum collapse time for this accuracy of detection which 
corresponds to the free particle limit with t co \ ~ 3.9/is. 
In current experiments detection efficiencies and band- 
widths will have a significant effect on the information 
that can be gathered from the record. For trapped ions 
harmonic frequencies would be around an order of magni- 
tude larger and state reduction times would then also be 
around an order of magnitude longer. However the cavity 
finesse used here for such experiments may be more dif- 
ficult to achieve at typical frequencies of ion transitions 
and the size of the cavity will be limited by the ion trap 
electrodes. What is dramatic about this time is that it 
is so short, single cold atoms have been observed close to 
the centre of the cavity mode in the experiment at Cal- 
tcch for times of the order of hundreds of microseconds. 
This is another confirmation of the extent to which this 
experiment operates near the Standard Quantum Limit 
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of position measurement ( |37j] and references therein) 
and an indication of the importance of continuous quan- 
tum measurement theory to its interpretation. 

Note that it will be difficult to achieve very small val- 
ues of r for which the measurement dynamics dominate 
over the oscillation in the well. This is the result of the 
two main constraints on the applicability of the model. 
Firstly it is necessary that the harmonic potential confine 
the atom to well within a wavelength, thus justifyin g th e 
Lamb-Dicke approximation in the master equation (2.1). 
this requires that the recoil frequency for the cavity tran- 
sition cj rec = M|/2m be much smaller than the harmonic 
oscillation. Moreover attaining the dispersive regime re- 
quires that the saturation parameter s = g 2 n/A 2 be 
much smaller than one. It is possible to express r in terms 
of these quantities: r — (uj/uj ICC )(l/8s)(uj/T cav ) where 
r^v ~_9q /k is the cavity mediated spontaneous emission 



rate 38 1 . Assuming that values of lj and s are chosen 
to satisfy a particular level of approximation, the rate 
at which the conditioned states become pure is only in- 
creased by changing the cavity parameters through in- 
creasing r cav relative to the oscillation frequency. More- 
over, since both the first two factors in this expression 
for r must be large for the SME (^T]) to correspond to 
the cavity system, extremely large values of r cav are nec- 
essary to achieve rapid measurement of the system state. 



III. CLASSICAL ANALOGUE 

If we are to interpret the conditioned state as the best 
description of the observer's knowledge of the quantum 
mechanical state given the results of a series of measure- 
ments, we would expect a similarity between these equa- 
tions and classical Bayesian state observation. The anal- 
ogy between the SME (2.1) and Kalman filtering for a 
classical position measurement was discussed in B but 
only the equations for the position probability distribu- 
tions were considered. Here we formulate the continuous 
time position measurement state observer for a classical 
harmonic oscillator and find that there is a close analogy 
between the SME and the classical theory for all mo- 
ments of the conditioned probability distribution as long 
as we restrict ourselves to Gaussian states and allow for 
noisy driving of the classical oscillator. 

The problem of noisy, classical, continuous time posi- 
tion measurement of a harmonic oscillator can be formu- 
lated 



dxc 
~~dt 

dpc 
dt 

dQ c 
dt 

E[e(t)e(t')\ 
E[e(t)r,(t')] 



-LOX C - 



ax c 



2^' 



E[rj(t)r)(t')}=5(t-t'), 
fS(t-t') 



(3.1a) 
(3.1b) 

(3.1c) 

(3.1d) 
(3.1e) 



We have used the same scaling of the variables as in the 
quantum problem. We imagine that as well as having 
an imperfect measurement of the system the oscillator is 
subject to a white noise force. There may be some cor- 
relation between the oscillator (plant or process) noise 
and the measurement noise and so e and r\ are correlated 
Wiener processes. As in the quantum mechanical case 
the limit of small r is the limit of good position measure- 
ment. 

The continuous time theory of Kalman filtering then 
provides the best estimate of the system state x,p and 
the second-order moments of the posterior probability 
distribution P(xc,Pc) 1 39 1 , 



dx c = ujpcdt + a\ — V xx dW, 



(3.2a) 



/ 2u 2u>\ 
d Pc = -tuxedt + aJ—V xp + fJ — dW, (3.2b) 



1 dVr- 



2V XV - -a'K 



2x^2 



u dt xp ' xx > 

oj dt \ \/rs i t 



2/ 2 ^ 2 2 



- -a'V„ 



Xp ' 



(3.2c) 



(3. 2d) 



j— = Vpp - 1 + —7= V xx a V xx V xp . (3.2e) 

lo dt \ Jrs it 



Where dW is an independent Wiener increment with 
dW 2 = dt proportional to the innovation process dQ~ax. 
Note that the circumflex employed here indicates that 
the quantity is an estimate of the classical variable and 
not that it is a quantum operator. The moments of 
the posterior probability distribution have been given 
the same notation as the moments of the conditioned 
quantum mechanical state. If we make the identifications 
s = r, a = cos ((>,/ = — sin <p then this system of equations 
is formally identical to the system which determines the 
evolut ion of the quantum mechanical conditioned state 
Eqns (Q. 

We see that there is a classical model of noisy po- 
sition measurement for which the equations of motion 
for the posterior probability distribution of the classical 
state given by the Kalman filter reproduce the stochas- 
tic master equation. What is specifically quantum me- 
chanical in the SME is that we cannot, even in princi- 
ple, specify the process noise and the measurement noise 
separately. Classically one could imagine isolating the 
system sufficiently that s is as large as we like. With 
s — ► 00 there would be no momentum disturbance on the 
atom and after a sufficiently long observation time the 
state of the system would be determined exactly so that 
Vl* = 0. Clearly this does not correspond to a quan- 
tum state. However the quantum theory of the prob- 
lem guarantees that any coupling to the system which 
gives position information about the state of the system 
must also disturb the momentum. This momentum dis- 
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turbance must be sufficient that the conditioned state 
always obeys the Heisenberg uncertainty principle. Thus 
only some classically possible models of position mea- 
surement are allowed by quantum mechanics. For a given 
level of measurement noise, the process noise must be at 
least sufficient to ensure that the observer can never infer 
probability distributions for the position and momentum 
which do not constitute valid quantum states. This is 
a measurement-disturbance uncertainty relation for 
continuous measurement; reducing the noise in the mea- 
surement must increase the noise in the evolution. This 
backaction noise does not however behave like classical 
process noise, its properties are entirely determined by 
the measurement. If we vary the basis of the measure- 
ment on the meter — vary <f> — then we are able to alter 
the correlation between the measurement noise and the 
apparent process noise in the classical model. 

We have found here that the symmetric moments of 
the conditioned state always obey a system of equa- 
tions which also describes Kalman filtration for a classical 
problem, but that altering the specific quantum measure- 
ment alters both the observation and the process noise 
of the relevant classical problem. Since the symmetric 
moments are the moments of the Wigner function, for 
example |27| , the conditioned Wigner function can be in- 
terpreted as the direct analogue of a classical posterior 
probability distribution for the system. This relationship 
is not going to be as straightforward for more compli- 
cated quantum systems where the conditioned Wigner 
function can be negative and cannot be interpreted as a 
probability distribution. 



IV. DEPENDENCE ON INITIAL STATE 

We have established above that the conditioned state 
of the system can in priniciple become pure after a finite 
observation time. What might not be clear is that the 
state which results is uniquely determined by the mea- 
surement record. In an ideal projective measurement 
the probabilities of obtaining the various measurement 
results depend upon the initial state. However, once a 
result has been obtained the conditional state after the 
measurement depends on that result alone. In a simi- 
lar way, for continuous measurements we find that, while 
the initial state affects the probability of obtaining par- 
ticular measurement records, the conditional state fol- 
lowing the measurement is determined by the particu- 
lar measurement record which was obtained, provided 
that record sufficiently long. In particular if the initial 
state of the system is very poorly known to an observer 
then we might hope that there is effectively a maxi- 
mum likelihood estimate of the system state which de- 
pends only on the measurement current and which con- 
verges to the actual system state within r co i. If prop- 
agating the stochastic master equation with the actual 
initial state of the system provides the a posteriori es- 



timate of the system state given the measurement re- 
sults, then such a maximum likelihood estimate would 
result from propagating an initial state with v ery large 
position and momentum variances in the SME (2.1), say 



V xx = V pp = Vo — * oo, V xp = 0, which gives a nearly flat 
prior probability distribution for the initial state. In this 
section we demonstrate that such a strategy does indeed 
work. Thus the purity of the conditioned density matrix 
indicates, as would be hoped, that there is only one pure 
state of the system at time t > t co \ which is consistent 
with the known sequence of measurement results. 

Suppose that two observers, Alice and Bob say, pos- 
tulate different initial states of the system p A (0), p B (0), 
which we will continue to assume to be Gaussian. For 
example, Alice may have more information than Bob 
about the system in which case Bob would start with 
a more mixed initial density matrix reflecting his initial 
lack of knowledge. Given that they both recieve the same 
measurement record I(t) and propagate their conditioned 
states according to the SME ( |2.l| ), it should be the case 
that at some time t' , around r co i, Alice and Bob agree on 
the system state, so that p A (t') = p B (t'). From the pre- 
vious section we know that after the time t co \ the second 
order moments of both conditioned states will be equal 
V A = V B = V?? and so we focus on the equations for 
the first order moments of the conditioned state for each 
observer 



dx A 
dp A 
dx B 

dp B = - 
dQ = x A dt 



ojp A dt + J —V x A dW 7 



cjx A dt+ \l —V xp dW, 



u P »dt + \I^V*dW B 



ujx a dt + x l—V B dW B , 



dW = x a dt 



dW t 



(4.1a) 
(4.1b) 
(4.1c) 
(4.1d) 
(4.1e) 



In this section we will omit the tildes used above to in- 
dicate that we have scaled the position and momentum 
to the natural units for the harmonic oscillator. The 
stochastic increment dQ is the infinitesimal increment of 
the measured current which both Alice and Bob have ac- 
cess to. We can express the incremen t dW B in terms of 
the other quantities dW B = dW — ^j2uj/r(x B — x A )dt, 
and find stochastic differential equations for the differ- 



ences between the means e x 



X , G^- 



P 



de x 



de r 



2uj. 
r 



ujtpdt V xx e x dt 



iy B -v A )dw, 



-u[l + -V£ I e x dt 



■ - 

(4.2a) 
(4.2b) 
(4.2c) 
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B 

xp 



v A ) dw. 



(4.2d) 



The deterministic part of this system of equations de- 
scribes a damped harmonic oscillation for e x , e p where 
the damping and oscillation rates depend on the second- 
order moments of Bob's conditioned state. The damping 
in these equations is not present in an analogous equa- 



tor a free particle. The 
btained from the continu- 



tion for e x given by Mabuchi [ 3 
equation adopted in |37| is o 
ous limit of the repeated position measurement model 
of Caves and Milburn and does not contain a noise 
term in the stochastic differential equation for x. In fact 
the omission of this term is in error and if the contin- 
uous limit of the repeated position measurement model 
of Caves and Milburn is taken correctly then the noisy 
contribution to de T , which we obtained from the SME 



(2.1), is in fact present. It is the damping which results 



from this term which leads to all observers agreeing about 
the conditioned state after a sufficiently l ong observation 
time. Note that after the time t co \ Eqs (4.2) are in fact 
a system of ordinary differential equations since at that 
point the covariances of the two conditioned states are 
equal. The differences in the means then damp to the 
steady state values e x — e p = indicating that Alice 
and Bob do eventually agree about the state of the sys- 
tem regardless of their initial states. Thus we have shown 
that the conditioned state eventually depends only on the 
measurement record but not the time-scale over which 
this occurs. 

It is straightforward to use the Ito chain rule |2(| to 
find differential equations for the expectation values of 
the covariance matrix for the difference between the con- 
ditioned means of Alice and Bob 



= 2E[e x e p ] 

uj ax 

+ - (V B 

1 V 



-V*E[4] 
r 

v A ) 2 



™ (E[e 2 p ]) = -2E[e x e p ] - ^V x B p E[e x e p ] 



1 d 

lo dt 



+± (v B -v A y 

r rip "xpj 
(E[e x e p ]) = E[e 2 p ] - E[&] 



(4.3a) 



(4.3b) 



2 V B E\el] - -V B E[e x e p ] 



xp 



(V B -V A X )(V B -V A ). (4.3c) 



Although we have taken the expectation value for this 
system of equations, the noise terms all become zero af- 
ter t co i so these ordinary differential equations eventu- 
ally describe the whole dynamics. We are now interested 
in the time-scale over which the determinant of this co- 
variance matrix becomes zero. Unfortunately the time 
dependence of the conditioned state variances prevents 
a closed form solution of this system of equations and 



we have not found a matrix Riccati form for the overall 
system. It is however straightforward to investigate the 
problem numerically. 

We wish to show that all observers will agree about 
the conditioned state of the system in roughly t co \. To 
do this it is sufficent to show that an arbitary observer 
will agree with some preferred observer in that time. For 
this reason we will assume that Alice has sufficient in- 
formation to describe the state of the system as a pure 
state and that she has access to sufficient earlier mea- 
surement records that V£{t) = V[*. An experimenter is 
going to be in Bob's position of not having precise knowl- 
edge of or control over the preparation of the state. We 
would expect that if Bob makes an accurate assessment 
of his initial knowledge of the state then E[e x ] is at first 
of the same order as V xx or smaller. Since V xx — V xx 
may be la rge initially, the stochastic term in the equa- 
tions ( |4.2j ) for Bob's errors dominates for very short times 
and essentially determines a random initial condition for 
e x such that E[e x ] < V xx . This reflects the fact that 
for short times the measurement record is dominated 
by noise. From this point on we find numerically that 
E[e x ] ~ V xx and that all the elements of the covariance 
matrix damp to zero within t co \. Thus we have found 
numerically that the conditioned state depends only on 
the measurement record after a time equal to the time 
over which the conditioned state becomes pure. Most im- 
portantly we found this to be the case even where all the 
variances were set to very large initial values, the largest 
we tried being V* = V* = E[e x ](Q) = E[e 2 p ](0) = 10 10 . 
Thus as we anticipated, Bob can make an accurate es- 
timate of the state within the collapse time even in the 
absence of accurate information about the initial state. 
In this case Bob's conditioned state corresponds to the 
maximum likelihood estimate discussed above. Interest- 
ingly, even with the pessimistic initial condition where 
_E[e^](0) is significantly larger than V B X (0) — this cor- 
responds to Bob overestimating the accuracy of his esti- 
mate of the particle's position and the difference between 
his estimate and the actual value being very large — the 
determinant of the covariance matrix damps on the char- 
acteristic timescale t co \ . This feature is more pronounced 
for larger values of V^(0). Again this is because for large 
initial values of V B , Bob's estimate of the system state is 
essentially dependent on the measurement current alone. 
These behaviours are demonstrated for the parameters 
used in the previous section in Figure (||). 

From this we can conclude that after the time r co i, any 
experimenter knows the state of the system regardless of 
how the system is initially prepared and of how much con- 
trol the experimenter has over this process. In the next 
section when we introduce detection efficiency and ther- 
mal couplings, these will just modify the second-order 
moments of the conditioned state such that the steady 
conditioned state is no longer pure. Since the equations 
for the first order moments will be unmodified all that is 
necessary to reproduce the results of this section in this 
more general case is to use the collapse time t co \ which is 







appropriate for the new system. Although this will mean 
that the experimenter is left with broader position and 
momentum probability distributions, it will not mean 
that different experimenters disagree about the means 
of these distributions. So regardless of the detection effi- 
ciency the conditioned state is eventually uniquely deter- 
mined by the measurement record and can be regarded 
as known by the experimenter. 



V. THERMAL AND DETECTION EFFICIENCY 
EFFECTS 

In order to discuss the effects of detection efficiency 
and other uncontrolled coupling to a bath in the model 
we will add an extra momentum diffusion term to the 
SME to obtain for 6 = 



dp c = —i[H sysl p c ]dt + 2otD[x]p c dt + 2j3V[x]p c dt 



2aH[x]p c dW. 



(5.1) 



This simple modification to the master equation is in- 
tended to model several possible imperfections in a real 
experiment. One contribution to /3 is the effect of detec- 
tion efficiency || so that (3 is at least a(l — r))/rj where 
the overall detection efficiency is r\. The effect of cavity 
loss through the unmonitored mirror in a cavity QED ex- 
periment is an effective detection inefficiency. Where the 
loss rates out of the two mirrors are Ki, /C2, and if only the 
light passing through the second mirror is detected, then 
we get (3 > otKi/pi where k = n\ + K2- The best situation 
is if the mirror in front of the detection apparatus has 
significantly higher transmission than the mirror used to 
drive the cavity. Scattering losses of the mirror will also 
be an effective detection inefficiency but these are typi- 
cally much smaller than transmission losses. In experi- 
ments with single atoms there will also be heating due 
to spontaneous emission into free space which will lead 
to a contribution to (3 equal to g^nT /4A 2 where T is the 
free space decay rate of the excited state of the atomic 
transition. Other than the restriction to the Lamb-Dickc 
regime this is the largest correction to the adiabatically 
eliminated master equation given above when a moderate 
detuning from the atomic transition is employed as in |Q. 
The spontaneous emission contribution to the heating is 
also proportional to the measurement coupling a such 
that j3 s = aTn/AgQ and to minimize the effect of spon- 
taneous emission we must use cavities with the largest 
possible value of g^/n. Recall that if this rate is large 
the signal to noise of the position measurement also im- 
proves. Contributions to (3 in other systems, for example 
the interferometric detection of the position of a moving 
mirror |l5|], will also come from any coupling of the os- 
cillator to a thermal bath. We found that the standard 
quantum Brownian motion master equation p?| ] led to 
steady conditioned states that did not obey the Hciscn- 
berg uncertainty principle for small values of r. This a 
result of the non-Lindblad terms in this master equation. 



The master equation we adopt here solves this problem 
by considering only coupling to very high temperature 
thermal baths for which the thermal contribution to (3 
is 7 >a , where 7 is the coupling rate to the thermal reser- 
voir and T is the temperature of the bath. This should 
be an adequate description of heating in the experiment 
as long as the bath to which the system is coupled is 
of sufficiently high temperature that only the diffusive 
evolution is significant for the timescales of inte rest . 

There is now no SSE equivalent to this SME (5.1) and 
pure states become mixed during the evolution. How- 
ever the SME continues to preserve Gaussian states so 
the previous calculation for the evolution of state can be 
straightforwardly modified. The second-order moments 
of the conditioned state still reach a steady state and we 
can easily find an expression for the steady conditioned 
state phase space area A is given by 



A ss = y/T+p/a = q > 1/0?, 
S (pf ) = 1 - 1/q > 1 - y/ri, 



(5.2a) 
(5.2b) 



S(PT) = ^ HQ + 1) - ^ 1* - 1) - ^ 2. (5.2c) 

The linear and von Neumann entropies of the steady con- 
ditioned states are plotted for a range of detection ef- 
ficiencies in Fig. (|j|). Even though we have effectively 
coupled the oscillator to an infinite temperature bath 
the conditioned steady states have in some sense a fi- 
nite temperature but stochastically varying mean values 
of position and momentum. 

We will consider some limiting cases here for the steady 
state variances in this more general situation. Where the 
measurement is strong (r <C 1) the position variance is 
insensitive to these imperfections, V* x ~ \fr^/q, ~ 
g, Vp* ~ 2y/q 3 /^/r. If the dynamics dominate (r 3> 1) 
then the position and momentum variances have the 
same dependence on (3 as A, V xx — q, V X p — q 2 /r, VJ^ ~ 

q- 

Finally the whole time evolution of the second-order 
moments of the conditioned state can be determined by 
solving the matrix Riccati equation. Again we will just 
consider the position variance as a function of time where 
14* (0) = V pp {0) = V >1, V xp (0) = 



V xx (t) 



q 2 V (c 2 cosh 2bt + b 2 cos 2ct) 
+q (V«? + q 2 ) (csinh26t - 6sin2ct) 



q 2 (b 2 + c 2 ) + qV (c 3 sinh 2bt + b 3 sin 2ct) 




(c 2 sinh bt — b 2 sin 2 ct) 



c = 7!VV 1 + 7? + 1 = ^ 



(5.3a) 



(5.3b) 



(5.3c) 



and we have scaled time by the harmonic oscillation fre- 
quency lo. As before when 2bt 3> lthe system reaches 
its steady state value and the non-linearity of the terms 
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describing the conditioning of the system state mean 
that the time for this to occur is independent of the ini- 
tial state. Interestingly this time is in fact shorter than 
was required to purify the conditioned state in the case 
of ideal detection. This is essentially because the ex- 
tra noise means that past observations become irrelevant 
more quickly, not leaving enough time to determine a 
pure state completely. While the steady state is reached 
increasingly fast it corresponds to an increasingly high 
effective temperature. When r > 1 this time to reach 
the steady state is r s ~ 2r/quj — muj/qha — T co \/q since 
in this regime b ~ q/r. As was noted in the previous 
section the time for the conditioned state variances to 
reach their steady state is also the time that is necessary 
for different observers to agree about the system state. 
The time evolution of some of these quantities is plotted 
in Figure (||). 

VI. CONCLUSIONS 

In this paper we have established that for a simple 
class of systems quantum trajectory theories allow the 
determination of a unique post-measurement state that 
depends only on the measurement results over a finite 
time. We have discussed the effects of experimental im- 
perfections on this state determination and the analogy 
to classical state observation. The systems to which this 
analysis is applicable have the property that the condi- 
tioned density matrix is at all times Gaussian and its evo- 
lution is exactly that of the posterior probability distri- 
bution for an appropriate classical state observer. While 
we have considered only the case of position measurement 
the same treatment will be applicable to these other lin- 
ear systems. Clearly in more complicated systems, such 
as the resonant interaction of an atom with the single 
mode of an optical cavity, this will not always be the 
case and the conditioned Wigner function will sometimes 
be non-positive. We would however expect that the cen- 
tral result we have shown here, the purity of the con- 
ditioned state after sufficiently long continuous observa- 
tion and the dependence of this state on the initial state 
only through the measurement results, will still hold for 
these more interesting and complicated systems. This 
will however require numerical simulation of the stochas- 
tic master equation for such systems. Another feature 
that should generalize is the interpretation of the SME 
as a state observer — presumably an optimal one — for 
the quantum system. 
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FIG. 4. (a) The time evolution of the unitless second-order 
moments V xx ,V xp ,V pp under continuous imperfect position 
measurement as described in the text. Time is measured 
in units of the harmonic oscillator angular frequency while 
r = 20 and q = 5. The initial state is V xx = V pp = 20, V xp = 
corresponding to a thermal state of the oscillator and for 
comparison the evolution of the conditioned position vari- 
ance is plotted for the case of perfect detection. As dis- 
cussed in the text the conditioned second-order moments 
reach steady state more rapidly when the detection is im- 
perfect, (b) The mean-square error in Bob's estimate 
of the conditioned state mean position and Bob's condi- 
tioned state position variance for imperfect continuous po- 
sition measurement r = 20, q = 5. The initial values are 
E[e 2 x ] = E[e 2 p ] = 10 10 , E[e x e p ] = 0, V X B X = V B = 10 10 , V X B P = 0. 
Even where the observation is imperfect the conditioned state 
is eventually independent of the initial state of the system. 



FIG. 1. The time evolution of (a) the unitless second-order 
moments V xx , V xp , V pp and (b) the two entropies of the condi- 
tioned state of the harmonic oscillator under continuous posi- 
tion measurement as described in the text. Note the different 
time axes for the entropies. Time is measured in units of 
the harmonic oscillator angular frequency while r = 20. The 
initial state is V xx = V pp = 20, V xp = corresponding to a 
thermal state of the oscillator. 



FIG. 2. The time evolution of Bob's mean-squared error 
in position E[e 2 ], Bob's conditioned state variance V B X is 
also plotted for comparison. In (a) the initial values are 
E[e 2 x ] = E[e 2 p ] = 5,E[e x e p ] = 0, V* = V* = 200,V x B p = 
and as discussed in the text the mean squared error, al- 
though originally small, rapidly becomes of the same order 
as the conditioned state variance. In (b) the initial values are 
E[e 2 x ] = E[e 2 p ] = 10 10 , E[e x e p ) = 0, V x \ = V% = 10 10 , V B = 
and over a time roughly equal to r co i the mean squared error 
approaches zero and the conditioned position variance ap- 
proaches its steady state. Note the different time-scales in 
the two graphs. 



FIG. 3. The linear and von Neumann entropies of the 
steady conditioned state are plotted against effective detec- 
tion efficiency. Extraneous heating of the oscillator is also 
described by r\ as discussed in the text. 
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